Crystalline and quasicrystalline local order in a model polymer melt 
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We investigate by Molecular-Dynamics methods the local order in a melt of a fully-flexible polymer chains. 
Evidence is given of clusters with well-defined crystalline and quasi-crystalline order. 
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I. INTRODUCTION 

There is growing interest in the possible presence of lo- 
calized ordered structures in liquids and polymeric systems 
JH Here we present the preliminary report of a thorough 
numerical study of these issues in polymer melts. By care- 
ful analysis of both the positional and the angular features 
of the inherent structures [3J, evidence is given of clusters 
with well-defined crystalline and quasi-crystalline order, i.e. 
face-centered cubic (fee) , hexagonally close-packed (hep) 
and icosahedral (icos) order. 

The paper is organized as follows. In Sec. |II]the polymer 
model and the Molecular-Dynamics (MD) algorithms are pre- 
sented. The results are discussed in Sec|IIIl Finally, the main 
conclusions are summarized in Sec II VI 



II. METHODS 

A. Details of simulations 

A coarse-grained model of a melt of linear fully-flexible 
unentangled polymer chains with fixed bond lenght and M 
monomers each is used. The system has N — 500 monomers 
in all cases but AI — 3, where N = 501. Nonbonded 
monomers at a distance r belonging to the same or different 
chains interact via a truncated LJ potential: 



ULjir) = e 
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+ Ucu 



(1) 



where a* = 2i/V is the position of the potential minimum 
with depth e, and the value of the constant Ucut is chosen to 
ensure Uq^p{r) — at r > Vc = 2.5 cr. Bonded monomers 
are costrained to a distance b — 0.97 cr by using the RATTLE 
algorithm All quantities are in reduced units: lenght in 
units of cr, temperature in units of e/fc^ and time in units of 
cr^ /i/e where ji is the monomer mass. We set fi = ks = 1- 
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States with different values of the temperature T, the density 
p and the chain lenght M; each state is labelled by the multi- 
plet {M, p, T}. NPT and NTV ensembles have been used 
for equilibration runs, while NVE ensemble has been used 
for production runs for a given state point. NPT and NTV 
ensembles are studied by the extended system method intro- 
duced by Andersen [5] and Nose |6]. The numerical inte- 
gration of the augmented Hamiltonian is performed through 
the multiple time steps algorithm, reversible Reference Sys- 
tem Propagator Algotithm (r-RESPA), developed by Tucker- 
man et at [7]. In particular, the NPT and NTV operators is 
factorized using the Trotter theorem separating the short 
range and long range contributions of the potential, according 
to the Weeks-Chandler- Andersen (WCA) decomposition J^. 



B. Inherent configurations 

At non-zero temperature monomers vibrate around their 
equilibrium positions; such fast movements make it difficult 
to characterize the arrangement of monomers. In order to re- 
move vibrations one usually resorts to the so-called inherent 
structures (IS) by mapping the configurations of the simulated 
trajectory into the corresponding local minimum of the po- 
tential energy |3]. The conjugate-gradient method is used to 
minimize the configurational energy as a function of the 3A^ 
particles coordinates |10]. Henceforth, a physical quantity X 
will be denoted as X^^ if evaluated in terms of IS configura- 
tions. 



III. RESULTS AND DISCUSSION 
A. Radial distribution function 

First, the radial distribution function g{r) of the system is 
considered. In a polymer system the function g{r) can be nat- 
urally separated into the intrachain part gic{r) and the inter- 
chain part gex{r)- 

In Figure [T] the intrachain part of the radial distribution 
function gff{r) is shown for trimers (M — 3) at differ- 
ent temperatures. The less-resolved instantaneous g{r) is 
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FIG. 1: Intrachain part of the inherent radial distribution function 
glc {r) for trimers (M = 3) with density p = 0.984 at the indicated 
temperatures. Inset: corresponding instantaneous intrachain radial 
distribution function gic{r). 
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FIG. 2: Monomer arrangements corresponding to the three peaks of 
the intrachain part of the radial distribution function in Fig[T] Dis- 
tance r and bond-bond angle 6 refer to the dark monomers, a) folded 
conformation (ri « 1.12, 6i ~ 70°); b) partially folded conforma- 
tion (r2 ~ 1.69, 02 ~ 120°); c) linear conformation (ra ~ 1.94, 
03 180°). 



also plotted. Bonded monomers give a S-liks contribution at 
r = b = 0.97. For larger distances, three peaks at ri, r2 
and ra are apparent. Figure |2] shows the arrangements of the 
three bonded monomers resulting in the distances ri, r2 and 
r^. The first and highest peak of gf^f (r) at ri « 1.12 coiTe- 
sponds to folded trimers in which the two chain ends are at a 
distance ri = \/2 w 1.12, the position of the minimum of the 
Lennard-Jones potential. The central peak at r2 ~ 1.69 cor- 
responds to the partially-folded conformation with two con- 
secutive bonds forming an angle 9 w 120°. The third peak at 
ra w 2fe w 1.94 signals trimers with linear conformations. 

The conformations plotted in Fig|2] are observed also in 
longer chain (not shown). 

In the panel (a) of Figure [3] the interchain part of the ra- 
dial distribution function gl^{r) is shown for trimers at dif- 
ferent temperatures. The most interesting feature is the split 
of the second-neighbors peak in three sub-peaks (enlargement 
in panel (b)), a feature which is apparent only in IS config- 
urations and is missing in the instantaneous configurations 



FIG. 3: Panel (a): Interchain part of the inherent radial distribution 
function gff (r) for trimers M — 3 with density p — 0.984 at the in- 
dicated temperatures. Inset: corresponding instantaneous interchain 
part of the radial distribution function gex{r). Panel (b): Enlarge- 
ment of the second-neighbors peak of gii{r). Note that the peak 
is resolved only by considering IS configurations (compare with the 
inset of panel a)). 
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FIG. 4: Monomer arrangements corresponding to the three subpeaks 
of the second-neighbors peak of gix(r). The distance is measured 
between the dark monomers. Bonds (black straight lines) are not 
shown if not relevant. Adjacent triangle (a') and face-sharing tetra- 
hedra (a") (r4 ~ r'l « 1.9); b): three linear monomers, two bonded 
(rs ~ 2.05); c): three linear nonbonded monomers (rg ~ 2.2). 



(Figj3^, inset). The arrangements of the monomers coiTe- 
sponding to the three sub-peaks of the second-neighbors peak 
of gi^{r) are shown in Figure H) The two side sub-peaks 
at r4 Ri 1.9 and rg w 2.2, already observed in atomic sys- 
tems ifTlll . signal two monomers separated by other overlap- 
ping monomers and three monomers in a linear conformation, 
respectively. In both cases the arrangements involve only non- 
bonded monomers. The central sub-peak at « 2.05 is pe- 
culiar of the polymer system and marks a linear arrangement 
of three monomers, two of them being bonded. For longer 
chains, the two side peaks vanish while the central peak in- 
creases slightly (not shown). 
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FIG. 5: Inherent distribution of the bond-bond angle P (cos 9) for 
trimers with density p = 0.984 at the indicated temperatures. Inset: 
corresponding instantaneous distribution P(cos 9). 



B. Bond-bond angle distribution 



0,6 

0,2 

0,6 



(icos) 
















1 








Further information on the local structure is gained by the 
distribution of the angle between adjacent bonds in a chain, 
P{cosd). In Figure|5]the distribution P^^{cosd) for trimers 
is shown. Again, IS configurations yield better resolution with 
respect to the instantaneous ones (Fig|5] inset). The distribu- 
tion vanishes above cos 9 ss 0.5 since monomers cannot over- 
lap. It is easily seen that the peaks of P^'^(cos 9) at 9i w 70°, 
62 ~ 120° and 6*3 w 180° originate from the dominant con- 
formations evidenced by the intrachain part of the radial dis- 
tribution function (Figure |2]l. 



C. Orientational order 

The radial pair distribution function g{r) provides only 
orientationally-averaged information on the arrangements of 
the monomers surrounding a given one. To gain insight into 
their orientational distribution, we resort to the rotational in- 
variant measures Qi jl2]. 

We consider the spherical polar angles 9{r) and 4>{r) of 
the vector r joining a central monomer with one of the rif, 
neighbors within a cutoff distance r* « 1.25a*. A spherical 
harmonic Yi„i[d{r),(p{r)] is associated to these angles. The 
quantity Qzm is defined as: 



rib 



(2) 



FIG. 6: Orientational order parameters Qi for 13-particles clusters 
with pure icosahedral (icos), face-centered cubic (fee) and hexago- 
nally close-packed (hep) order. 
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FIG. 7: Orientational order parameters Q'^ for IS configurations of 
trimers (M = 3) and decamers (M = 10) at density p = 0.984 and 
temperature T = 0.6. 



The orientational order parameter Qi is then defined as: 



(3) 



In the absence of orientational order Qi — 0, V/ 7^ 0. The 
Qi values for different kinds of orientational order are easily 
computed. In Figure |6] several order parameters are shown 



for three clusters of 13 particles with pure icosahedral (icos), 
face centered cubic (fee) and hexagonally close-packed (hep) 
order. Non-zero values appear for I ^ A in the fee and hep 
clusters while in the icosahedral cluster non-zero values occur 
only for I — 6 and I — 10. 

The orientational order parameters for I — 

2, 4, 6, 8, 10 are shown in Figure [T] for IS configurations of 
trimers Af = 3 (upper panel) and decamers M — 10 (lower 
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panel). The presence of icosahedral orientational order is sig- 
naled by Qio 7^ 0. However, the non-zero values of Qi for 
Z = 4, 8, which vanish under icosahedral symmetry, indicate 
the presence of partial fee and hep orientational order too. 

It is easily seen that the positional and the angular features 
evidenced by the analysis of the radial distribution function 
and the bond-bond angle distribution in Sections IIII Al and 
HUB I are further evidence of the different kinds of orienta- 
tional order shown by the orientational order parameters. Fi- 
nally, one notes that the order parameters Qi are coincident 
for trimers and decamers, thus suggesting that the increasing 
connectivity does not affect the local order in a marked way. 



IV. CONCLUSIONS 



By analysing the positional and the angular features of the 
IS radial distribution function and bond-bond angle distribu- 
tion as well as the local orientational order parameters, we ev- 
idenced the presence of clusters with crystalline (hep and fee) 
and quasi-crystalline order in a melt of fully-flexible polymer 
chains. 
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